Anisotropy of the Optimally— Doped Iron Pnictide Superconductor 

Ba(Feo.926COo.074)2As2 
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Anisotropies of electrical resistivity, upper critical field, London penetration depth and critical 
currents have been measured in single crystals of the optimally doped iron pnictide superconductor 
Ba(Fei_i,Coi,)2As2, 2;=0.074 and Tc ~23 K. The normal state resistivity anisotropy was obtained by 
employing both the Montgomery technique and direct measurements on samples cut along principal 
crystallographic directions. The ratio 7p = pal pa is about 4±1 just above Tc and becomes half of that 
at room temperature. The anisotropy of the upper critical field, = IIc2,ab/ Hc2,c, as determined 
from specific heat measurements close to Tc, is in the range of 2.1 to 2.6, depending on the criterion 
used. A comparable low anisotropy of the London penetration depth, 7^ = Xc/Xab, was recorded 
from TDR measurements and found to persist deep into the superconducting state. An anisotropy 
of comparable magnitude was also found in the critical currents, 7^ = jc,ab/jc,c, as determined from 
both direct transport measurements (~1.5) and from the analysis of the magnetization data (~3). 
Overall, our results show that iron pnictide superconductors manifest anisotropies consistent with 
essentially three-dimensional intermetallic compound and bear little resemblance to cuprates. 



INTRODUCTION 

Discovery of the iron-arsenide family of high crit- 
ical temperature superconductors Q naturally raises 
the question about the relation between them and the 
cuprates These materials share many common fea- 
tures. Both materials are layered with electronically ac- 
tive Cu-0 and Fe-As layers alternating with buffer layers 
of different chemical composition. The 3d electronic or- 
bitals of copper and iron make the main contribution 
to the electronic bands close to the Fermi energy. As 
a result of this layered structure, the cuprates reveal 
highly anisotropic electronic properties, with the ratio 
of conductivities along and perpendicular to the con- 
ducting layer, 7p = pd Pa, varying from about 50 in 
YBasCuaOr H to above 10^ in BisSrsCaCuzOio [I] at 
optimal doping. This high anisotropy is a reflection of a 
two-dimensional Fermi surface, as found experimentally 
in semiclassical and magnetoquantum oscillations Q. 
High anisotropies are also found in the lower [3] and up- 
er critical fields [s*] and in the critical current density 
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For the iron arsenides, early band structure calcula- 
tions have also suggested a two-dimensional electronic 
structure fl^. High anisotropy of the resistivity, with 
7p ~100, was reported for the non-superconducting par- 
ent compounds BaFe2As2 fill and SrFe2As2 [i21j as well 
as for Co-doped BaFe2As2 [13|. Contrary to this, a rela- 
tively low anisotropy is found in the upper critical field of 
all studied iron arsenide compounds close to optimal dop- 
ing [ll [H, [H, [13, [H, [11, S [m, m, [H, m . The angular 
dependence of resistivity as a function of magnetic field 
also suggests a rather small anisotropy in NdFeAs(0,F) 
[2^. The anisotropy evolves with doping, showing a 



two-fold change between the underdoped {x <0.074) and 
overdoped {x >0.074) regions [2J]. This suggests that 
the anisotropy may be an important parameter to char- 
acterize superconductivity in the iron pnictides. 

The anisotropy of the upper critical field, 7/f = ^jpf^, 
and of the London penetration depth, 7;^ = Xc/Xab, 
are linked in the region of validity of the Ginzburg- 
Landau (GL) theory of phase transitions close to Tc, 
also for isotropic gap and in the dirty limit, ^ y/jp 
0, [23, m, [23I. Gross violation of these relations and 
very high values of Hc2 [lO] can indicate, for example, 
paramagnetic [3l| or some exotic [s^l mechanism of su- 
perconductivity suppression already in the very vicinity 
of Tc- This situation is realized in two-dimensional or- 
ganics for magnetic fields applied parallel to the super- 
conducting planes [s^, in the spin-triplet superconduc- 
tor Sr2Ru04 34 1 and in the heavy fermion CeCoIns [3^. 
In all of these cases the temperature interval of the va- 
lidity of GL theory is very small and the temperature 
dependence of jh at any sizable field does not follow GL 
predictions. 

In this work we have undertaken comprehensive char- 
acterization of the anisotropy of Ba(Fei_a;Coa;)2As2, 
x=0.074, with Tc w23 K. Our choice was motivated by 
the availability of high quality single crystals and very 
good reproducibility of their properties between differ- 
ent groups for all doping levels [23, ISfil, \3^. We have 
found that the anisotropies of the upper critical field, 
electrical resistivity, London penetration depth and crit- 
ical current, which were determined using different mea- 
surements, agree with each other and show values much 
lower than in the cuprates, indicating that the salient 
physics associated with superconductivity in these two 
families of compounds may be different. 
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EXPERIMENTAL 

Single crystals of Ba(Fei_a;Coa;)2As2 were grown from 
FeAs/CoAs flux from a starting load of metallic Ba, FeAs 
and CoAs mixed in the proportions 1:3.6:0.4, as described 
in detail elsewhere [24] . Crystals were thick platelets with 
sizes as big as 12x8x1 mm'^ and large faces correspond- 
ing to the tetragonal (001) plane. The cobalt content 
in the crystals was determined by wavelength dispersive 
X-ray electron probe microanalysis to be a;=0.074. The 
crystal quality of the samples was confirmed with X-ray 
Laue measurements on single crystals, which found reso- 
lution limited narrow peaks, see and [sl] for details. 

It was shown in previous study that correct determi- 
nation of the sample resistivity is not a simple problem 
for the iron arsenides 2j|. Due to softness of the ma- 



terials, their cutting and shaping into transport samples 
inevitably introduces cracks, which affect the effective 
geometric factors of the sample. This represents an es- 
pecially serious problem for measurements with current 
along the c axis. A strong tendency to exfoliate prevents 
the cutting of samples with c 3> a. As we will show 
later, partial cleaving by exfoliation is one of the most 
likely reasons for the unusu ally high anisotropy, as found 
m previous studies [lil,[i3,|l|. 

Samples for electrical resistivity and critical current 
measurements with current flow along the [100] a-axis in 
the tetragonal plane {pa and Jc,a) were cut into bars of 
5*0.12*0.025 mm3 {a x b x c), as described in Ref. 
Samples for the Montgomery technique measurements 
(see below) of the resistivity anisotropy ratio, 7p, as well 
as for electrical resistivity and critical current measure- 
ments with current flow along the tetragonal c axis (pc 
and Jc,c ) were cut into (0.2-0.5)*(0.2-0.5)*(0.1-0.5)mm3 
{a X b X c) bars. 

Contacts to the samples were made by attaching sil- 
ver wires with a silver alloy, resulting in an ultra low 
contact resistance (less than 100 p,n). Measurements of 
Pa and Jc^a were made in both standard 4-probe and 2- 
probe configurations and gave identical results, see (soj 
for details. Measurements of pc and Jc,c were made in 
the two-probe sample configuration. Even for the nomi- 
nal 2-probe measurements, a 4-probe scheme was used to 
measure the resistance down to the contact to the sam- 
ple, i.e. the sum of the actual sample resistance and 
contact resistance Rc was measured. Since Rg 3> Rc, 
this represents a minor correction. This can be directly 
seen at temperatures T < Tc, where Rg =0 and measured 
resistance represents Rc (see Fig. |4] and Fig. [5] below) . 

In 1961 Wasscher [i^, based on van der Pauw calcula- 
tions [41], found that the current distribution in a sample 
with dimensions of h and I2 along principal directions of 
the resistivity tensor, pi and p2, is equivalent to that 
of an isotropic sample with dimensions h and ^2\/ P2/ Pi 
4lj. This scaling transformation is used in the Mont- 
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FIG. 1: Inset: Arrangement of electrical contacts in the 4- 
probe Montgomery technique for determination of electrical 
resistivity anisotropy. The sample is cut into a rectangular 
prism, with the sides of its bases h and I2 oriented along the 
principal directions 1 and 2 of the conductivity tensor, a and 
c crystallographic directions in our case. Contacts (red lines) 
are put over the whole length of the sample in the third direc- 
tion, ^3 (II a). Two successive 4-probe measurements are made 
with current (blue arrows) flowing between contacts on one 
side of the sample and voltage (green arrows) measured on 
the other side of the sample, to determine resistance values 
Ri (current along li) and R2 (current along h). Main panel: 
Temperature dependence of resistances _Ri (small solid sym- 
bols) and R2 (large open symbols) for three different samples 
of Ba(Fei_^Co^)2As2, 2;=0.074. 



gomery technique to map measurements on samples of a 
known geometry and an unknown anisotropy onto those 
in isotropic samples, where this ratio can be calculated 
4^ . 43j . Two successive 4-probe resistance measurements 
are made using the contact configuration shown in the in- 
set in Fig. [TJ First, current flows along the li direction, 
in our case corresponding to the a-axis in the plane, pro- 
ducing a voltage drop on the opposite side of the sample 
to determine the resistance Ri = Vi/Ii. Second, the di- 
rection of the current is rotated by 90 degrees along the 
I2 direction (along the tetragonal c-axis in our experi- 
ment), and the resistance R2 is determined. The ratio of 
the measured resistances, R1/R2 is exponentially sensi- 
tive to the ratios of sample dimensions, Z1/Z2, and is used 
for determination of the resistivity anisotropy [42] . Since 
the whole idea of the Montgomery technique is based on 
a homogeneous current distribution in the sample, the 
structural integrity of the sample plays a crucial role for 
measurements of this type. 

Analysis of the measured anisotropy signal was per- 
formed assuming the precise position of the contacts at 
the corners of the sample and neglecting their size in 
the basal ac plane of the rectangular prism. The con- 
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tacts were extending along the sides perpendicular to the 
basal plane, so a thin slab approximation was used in the 
analysis. Since the actual size of the contacts is not neg- 
ligible as compared to the sample circumference in the 
ac plane, and their positions can deviate from ideal, this 
brings sizable errors into the anisotropy measurement. 
By reproducibility of the results on three measured sam- 
ples, we estimate the error to be on the order of ± 50% 
for the anisotropy ratio. 

Transport critical current measurements were per- 
formed by measuring IIV characteristics at fixed tem- 
peratures, see [3^ for details. Samples for Jc.a measure- 
ments were mounted with GE-varnish on an insulating 
oxide substrate to serve as a heat sink and contacted with 
either multiple silver wires or silver foil to create large 
area contacts providing an additional heat sink. Sam- 
ples for Jc,c measurements were sandwiched between two 
silver foils and thermally grounded to large silver heat 
sinks. Currents up to 2 A were generated in a commercial 
PPMS measuring system from Quantum Design. Mag- 
netic critical current densities were estimated by using 
the Bean model 4J| from the magnetization measure- 
ments performed in a Quantum Design SQUID magne- 
tometer and independently from the profiles of the mag- 
netic induction, measured by utilizing a magneto-optical 
method [ist . 

Measurements of the penetration depth were per- 
formed by using a tunnel diode resonator (TDR) tech- 
nique [4^ . This technique is capable of resolving changes 
in the penetration depth of about 1 A. Details of the mea- 
surement technique are described elsewhere 47[ . In brief, 
a properly biased tunnel diode compensates for losses in 
a tank circuit, so it is self-resonating at a frequency of 
/ = l/{2n\/LC) ~ 14 MHz. A sample is inserted into 
the coil on a sapphire rod. The change of effective induc- 
tance causes a change in the resonant frequency. This 
frequency shift is proportional to the dynamic magnetic 
susceptibility of the sample, x- Knowing the geometrical 
calibration factors of the circuit, we obtain X{T,H) as 
described in Ref. 47|. A major advantage of this tech- 
nique is a very small AC excitation field amplitude 20 
mOe), much lower than Hd ^ 100 Oe. This means that 
the TDR technique only probes, but does not disturb 
the superconducting state. Other advantages are high 
stability and excellent temperature resolution 1 mK). 

The upper critical field Hc2 was determined from the 
onset of the superconducting transition in the TDR mea- 
surements as well as from measurements of the specific 
heat, taken in a Quantum Design PPMS. The same sam- 
ple was used in the specific heat, DC magnetization mea- 
surements of the critical currents and in the magneto- 
optical imaging. In the configuration with H \\ ab, the 
sample was aligned parallel with the field using a pre- 
cisely cut sapphire substrate. Two samples were studied 
and they have shown similar properties. 

Band structure calculations have been done within the 



full potential linearized augmented plane wave (FLAPW) 
approach (48| within the local density approximation 
(LDA) [4§|. The mesh of 31x31x31 fc-points was used 
for the Brillouin zone integration and Fermi surface 
plot. We have used experimental lattice constants for 
the BaFeaAsz [H^ and Ba(Feo.926Coo.o74)2As2 0. The 
Fermi velocities were calculated using the Bolz-Trap 51 1 
package. 



RESULTS 
Anisotropic resistivity measurements 

In Fig. [T] we show the temperature dependence of the 
experimentally measured resistances Ri (current along 
the plane, solid symbols) and R2 (current perpendicu- 
lar to the plane, open symbols) for 3 different samples. 
At room temperature, Ri is higher than R2 in all the 
samples, as expected for samples of larger in the plane di- 
mension and small anisotropy. The resistivity anisotropy, 
7p = Pel Pa, was deduced from the Montgomery proce- 
dure and is shown in Fig. [51 There is an overall general 
agreement in the temperature dependence of both resis- 
tances for the three samples measured, however, with 
variation of R1/R2 due to its exponential dependence 
on the ratio of the sample dimensions. The anisotropy 
7p varies between 2.2 and 3.5 at room temperature and 
increases on cooling approximately by a factor of two, 
reaching (3 to 5) at Tc- Resistance jumps seen in two 
samples most likely indicate that the samples undergo 
partial cracking on cooling. Since the Montgomery tech- 
nique heavily relies on the idealized current distribution 
in the sample, it is impossible to ascertain that the ratio 
is determined correctly at temperatures below which the 
crack formation happens. However, at room temperature 
and at temperatures down to the appearance of cracks, 
the data seems to be quite reliable and well reproducible. 
Of the three samples, the most reliable measurement was 
done on sample #1 (black curves), where R\ and R2 are 
comparable over the entire temperature range, implying 
good compliance with requirements of Montgomery anal- 
ysis. Cracks affect both i?i and i?2, suggesting current 
redistribution in the sample, although resistance values 
do not change dramatically. 



Artifact of high anisotropy and different 
temperature dependence 

In Fig. [3] we show the Ri and R2 data from yet another 
sample, measured in the Montgomery configuration. At 
room temperature this sample shows R2 ~ 7 * de- 
spite comparable sample dimensions, suggestive of no- 
tably higher anisotropy, as can be seen in the inset of 
Fig. [21 Interestingly enough, Ri and R2 individually re- 
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FIG. 2: Temperature dependence of the ratio of resistivities 
pel Pa as determined from the Montgomery technique. For 
comparison we show the temperature dependence of the ratio 
as determined from independent pa and pc measurements (red 
line). 




FIG. 3: Left panel: Temperature dependence of the raw re- 
sistance data _Ri and K2 (main panel) and of the resistivity 
anisotropy ratio pc/pa (inset) for sample 4. Note very differ- 
ent temperature dependence of pa/ pc as compared to samples 
1-3, in Fig. (2] despite a very similar temperature dependence 
of individual values of R\ and -R2. Right panel: photograph of 
the sample after two month storage. On degradation, sample 
revealed clear crack perpendicular to the c-axis, which was 
not noticeable after sample preparation. 



veal a similar temperature dependence to those found in 
the other 3 samples, however, the ratio of the measured 
resistivities, as determined from the analysis, gives a very 
different temperature dependence as compared to sam- 
ples 1-3. As we will show later, the anisotropy ratio de- 
termined from direct resistivity measurements shows the 
same temperature dependence as the ratio determined 
from the Montgomery technique measurements on sam- 
ples 1-3. Based on this qualitative difference we conclude 
that the data for this sample is not representative. 

We should point, however, that even with such no- 
tably different measured Ri and i?2, the anisotropy ra- 
tio increases by a factor of 5 or so, but notably distorts 
the temperature dependence. This stresses the need for 
direct resistivity measurements with the current along 
principal directions of the conductivity tensor. 

To get insight into the possible reason for this behav- 
ior, we inspected the sample with unusual anisotropy ra- 
tio. The right panel of Fig. [3] shows the sample after two 
months storage. Due to sample degradation, the cracks, 
which were not noticeable originally, developed and be- 
came noticeable. This particular sample almost split into 
two pieces. 



Temperature dependence of in-plane and inter-plane 
resistivity 

Since Montgomery resistivity measurements allow for 
the correct determination of the ratio of resistivities, but 
not their individual temperature dependences, we need 
to perform direct measurements of at least one of the 
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FIG. 4: Temperature dependence of pa measured on 4 dif- 
ferent samples from 3 different batches. Left panel- actual 
measurements, right panel- the data normalized to the room 
temperature values. 

components. To check the consistency of the obtained 
results, we have independently measured the resistivity 
for both current flow directions. 

The data for Pa{T) are shown in Fig.[4l The data are 
taken on samples from 3 different batches. Note that the 
scatter of resistivity values at room temperature is far 
above the error bar in the geometric factor determina- 
tion. On the contrary, as shown in the right panel, the 
shape of the temperature dependence remains the same, 
with Pa{T) varying roughly linear in T over the entire 
temperature range between room temperature and T,.. 
The data can be actually matched even better allowing 
for slight variation of the residual resistivity po- This is 
consistent with the results of our study using different 
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FIG. 5: Temperature dependence of pc measured on 3 sam- 
ples. Left panel- actual measurements, right panel- the data 
normalized to the room temperature values. 



sample cutting/contact making technique [24 1. 

This similarity of the shape of the temperature depen- 
dence together with a large difference in the inferred resis- 
tivity values clearly shows that the problem comes from 
incorrect determination of the geometric factor. Two 
types of defects can be invoked to explain this discrep- 
ancy. Inclusions of FeAs flux usually happen between 
layered crystallites. Cracks tend to exfoliate layers. In 
both cases the effect will be more visible for c-axis trans- 
port, but will affect the Pa measurements as well. How- 
ever, since FeAs is a metal, it is hard to imagine that 
its own contribution to the conductivity will not affect 
the shape of the resistance temperature dependence. No 
evidence for FeAs flux inclusions is found also in x-ray 
studies [13. 

The temperature dependence of pc taken on three dif- 
ferent samples is shown in Fig. [Sj Due to small sample 
size along c-axis we were only able to perform two probe 
resistivity measurement. As discussed in Methods sec- 
tion above, the contact resistance is negligibly small. The 
three curves are similar, which is better seen when the 
data are normalized by room temperature values, Fig.[5l 
right panel. They show an extended range of weak vari- 
ation in pc{T) down to approximately 100 K followed by 
a roughly linear decrease below. A factor of 4 difference 
in the pc value implies variation of effective sample cross- 
section due to cracks. 

To test the possible effect of cracks on Pc{T) measure- 
ments we performed two successive measurements of the 
resistivity and of the critical current on the same sam- 
ple. After an initial run with Pc(T) and Jc measure- 
ment, we perceived that pc was potentially artificially 
high. To test this the sample was cleaved into two pieces 
by pulling contact wires and applying small force parallel 
to the plane. A new contact was made on fresh cleaved 
surface, and measurements were repeated. The resultant 
Pc{T) and V — J curves are compared in Fig. [S] As is 
clear from the figure, cleaving decreases resistivity of the 
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FIG. 6: Left panel Temperature dependence of pc measured 
on the same sample after initial preparation and after cleaving 
the part of the sample and remaking contacts. Right panel 
shows V-J curves for the same sample at 21.4 K, revealing 
increased critical current density after sample cleaving. 



sample (beyond changing only geometric factor) and in- 
creases critical current density. This is consistent with 
the sample initially having a crack perpendicular to the 
c-axis, giving rise to the high pc value. In addition, this 
clearly shows that samples with the lowest measured re- 
sistivity Pc should be used for evaluation of resistivity 
anisotropy. 

For Pa, the effect of cracks due to exfoliation is more 
complicated. Here cracks not only affect measured re- 
sistivity but can disrupt connection between current and 
voltage contacts. As a result this can either increase (if 
current flows across the crack) or decrease (if voltage con- 
tacts are disconnected from current path) measured re- 
sistance values Because of this, we have excluded 
from consideration the extreme curves. As an additional 
criterion for the selection, we have used the data for the 
sample with the highest critical current density, which 
would obviously exclude samples with cracks. 

In Fig. [2] we show a comparison of the temperature de- 
pendence of the anisotropy ratio from independent mea- 
surements of Pa{T) and pc{T), and the ratio as measured 
by using the Montgomery technique. Taking into account 
that these data are taken on different samples in very 
different measurement conditions and a sizable amount 
of uncertainty in geometric factors, the agreement be- 
tween the two independent anisotropy determinations is 
remarkable. 

To evaluate error bar for the determined anisotropy 
ratio, we need to estimate a spread in resistivity values 
for Pa and pc- As we have shown above, simple averag- 
ing of Pc is not meaningful. As extreme case scenarios 
we take minimum (maximum) resistivity pa and pc and 
come to the anisotropy pd Pa between 1.2 and 10 at room 
temperature. This covers all the range of Montgomery 
technique determinations, and is notably lower than the 
number reported in previous study [l3| . 
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FIG. 7: Temperature dependence of specific heat, represented 
on a C /T vs T plot, zooming into the range close to Tc for 
magnetic fields H \\ c (left panel) and HLc (right panel). 



Anisotropy of the upper critical field 
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FIG. 8: Temperature dependence of the difference of specific 
heat, C/T-C/T{H \\ c, 6 T), for H \\ c (left panel) and H_Lc 
(right panel). This subtraction removes non-superconducting 
background to the specific heat for T > Tci&T, H \\ c). The Tc 
was determined from the position of the maximum in the data 
and from the position of an onset, determined as a crossing 
point of linear extrapolations of the rising part of the C/T 
jump and of the constant data above Tc, as shown. 



In order to compare the anisotropy of the electrical 
resistivity to the anisotropy of iJc2, we need to deter- 
mine Hc2{T) near Tc- Two previous studies performed 



resistive determination of Hc2 [23, l2J] on samples of 
Ba(Fei_a;Coa;)2As2 with similar Tc and composition. In 
Ref. [Ill x=0.1 according to the starting load, similar 
to our samples, which actually corresponds to 0.0 74 | 24| . 
Both found low ^h- Close to T^ jh is about 2 [22| or 
2.5 to 3.2 [23|. On cooling much below Tc it decreases to 
1.5 at ~ Tc/2 or 1.5 to 2 at 0.7Tc The slopes, 
dHc2/dT, were found to be 4.5 T/K and 2.4 T/K for 
H \\ ab and H \\ c, respectively [3a], however, these are 
strongly criterion dependent [2J|. This projects to a resis- 
tivity anisotropy 7p ^ ^fj of 4 to 9, which is in reasonable 
agreement with our findings. However, resistive determi- 
nation of Hc2 is strongly criterion-dependent and poten- 
tially subject to contributions of superconducting fluctu- 
ations and current percolation. This can notably com- 
plicate determination of Hc2 at the lowest fields, where 
changes are small. With this in mind, we have decided 
to use specific heat measurements as a bulk probe of the 
Hc2i revealing simultaneously sharp and easy to recog- 
nize feature of the superconducting transition. 

Fig. [7| shows the temperature dependence of the spe- 
cific heat obtained on a m=1.9 mg sample, used in 
magneto-optical and magnetization studies of the criti- 
cal currents, Fig. [7| This sample is of very regular shape 
which allowed for its precise alignment in magnetic field. 

The jump in the specific heat at Tc broadens with the 
increase of magnetic field for both directions of an ap- 
plied magnetic field (Fig. [7|). This makes the determi- 
nation of Tc{H) more ambiguous. To obtain the best 
resolved jump at the lowest fields, which are of main in- 
terest for our comparison, we subtract the C/T data for 
H{\\ c) — from all low field curves, as shown in Fig.[5| 
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FIG. 9: Temperature dependence of the upper critical fields 
as determined from the onset (triangles) and the maximum 
position (stars) in specific heat data and from TDR data (cir- 
cles). The data for H \\ ah is shown with solid symbols, for 
II c is shown as open symbols. 



This subtraction removes the large non-superconducting 
background to the specific heat and reveals a sharper 
C/T jump. To determine Tc{H) we have used linear fits 
of the rising portion of the C/T data, as shown in Fig. [5| 
Thus we have determined the onset of the specific heat 
anomaly, while the maximum position was used as yet an- 
other criterion for the Tc{H) determination. Depending 
on the criteria used, we get different slopes of the Hc2{T) 
at Tc. For H \\ ab [H \\ c) it ranges from 8.8 (3.4) T/K 
from the onset of the C/T jump to 5.3 (2.45) T/K from 
the position of the maximum. 
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FIG. 10: Temperature dependence of the TDR frequency 
shift in magnetic fields applied along the tetragonal c-axis (left 
panel, -ff || c) and perpendicular to it (right panel, H \\ a). 
The superconducting transition temperature was defined as 
a crossing point of linear extrapolations of the shift in the 
normal state and at the slope. 

Fig. [9] shows a comparison between bulk thermody- 
namic Hc2 determined by both criteria from the specific 
heat with Hc2{T) determined from TDR measurements. 
The TDR data used for the Hc2 determination are shown 
in Fig. 1101 In zero apphed external magnetic field, the fre- 
quency shift is representative of the variation with tem- 
perature of the London penetration depth. In constant 
magnetic field it represents the onset of the supercon- 
ducting shielding and is close to the onset of the resistive 
transition fs^. The slopes determined from the TDR 
measurements are 4.56 T/K for H \\ ah and 1.87 T/K for 
H \\ c. Thus by all the criteria used we have obtained 
between 2.1 and 2.6, which matches within the error 
bars to the electrical resistivity anisotropy determined 
from the Montgomery technique. 



Anisotropy of the London penetration depth 

The in-plane penetration depth, AAa6(T), was deter- 
mined by applying the rf magnetic field along the c-axis 
( Hr/llc). In this geometry, the screening currents flow 
within the a6-plane, so only the in-plane component of 
the penetration depth is probed. When H^/ -L c, the 
screening currents flow both in-plane and along the c- 
axis. Using the model described by Ref. [43], the c-axis 
component AAc was obtained. Because the TDR tech- 
nique measures changes in penetration depth with tem- 
perature (AA(T)), in order to calculate A(T) we need 
as a reference, the absolute value of penetration depth 
at a certain temperature. For the in-plane component, 
we have used the zero temperature value, Aa6(0) « 208 
nm, from Ref. [s^ . As we are not aware of any previous 
reports of Ac(0), for the c-axis component we make the 
assumption that near Tc, according to Ginzburg-Landau 
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FIG. 11: Variation of the London penetration depth, 5\ = 
A(r) — A(r — 0) with temperature for two orientations of 
AC field with respect to the sample. Inset shows tempera- 
ture dependence of the anisotropy of penetration depths, 7a, 
assuming 7a = jh a,t T — O.QTc and two selected values of 
7i/= 2.0 and 2.5. 

theory, both the upper critical field and the penetration 
depth anisotropics should be equal. Considering 7^ = 7a, 
we set at T=0.9Tc, Ac = Jh x Aaf, with two selected val- 
ues for "fH= 2.0 and 2.5. The resulting low anisotropy 
and its weak temperature dependence are consistent with 
the overall picture of small anisotropy in the pnictides 

(Fig.ni. 

Anisotropy of the critical current 

In Fig. [T^l we show isothermal current density- volt- 
age, J-V, dependences for currents flowing along a-axis 
(left panel) and along tetragonal c-axis (right panel). 
The resistivity of the sample used in the critical cur- 
rent measurements along the c-axis at room temperature 
was about 800 /Ltf2cm, a value consistent with minimal 
cracking. For both directions of the current flow we have 
determined the same resistive transition with the onset 
at 22.5 K, the midpoint at 22 K and the zero-resistance 
state at 21.6 K (±0.05 K in all cases). The VJ curves are 
linear in the normal state, implying that Joule heating is 
insignificant. The critical current was determined at the 
point of the sharpest voltage rise, for which the deriva- 
tives of the actual V- J curves were taken [s^ . This proce- 
dure is unambiguous when the sample is cooled somewhat 
below a temperature where its resistance becomes zero. 
In the inset of Fig. [13] we compare the critical current 
densities determined from the transport measurements. 
The data reveal only modest anisotropy (less than 1.5). 
However, the temperature range for the determination 
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FIG. 12: Isothermal dependences ol voltage drop on the sam- 
ple, V, vs current density, J, taken in the vicinity of the resis- 
tive transition to the superconducting state. Left panel, cur- 
rent along o-axis, right panel- current along tetragonal c-axis. 
For current densities higher than certain values the curves be- 
come linear. The critical current density was determined as 
the position of the maximum of the derivative of V{J). 
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FIG. 13: Superconducting critical currents as determined 
from transport measurements (stars in the inset) and from 
magnetization measurements using the Bean model (main 
panel and inset, solid circles). 

of the transport critical current density is limited to the 
immediate vicinity of Tc, and small differences in Tc can 
potentially affect this ratio. To determine the critical 
current anisotropy in a broader temperature range we 
turn to the magnetic measurements. 

We use the Bean model with field- independent, con- 
stant throughout the sample, supercurrent density js- 
Assuming that js is isotropic in the crystallographic a6— 
plane, three different current densities are possible de- 
pending on the orientation and the direction of motion of 
the Abrikosov vortices under the influence of the Lorentz 
force. 

Let us consider a rectangular prism - shaped crystal 



of dimensions, 2c < 2a < 2b. In the crystals of interest, 
the crystallographic a6— plane has the largest area and 
is parallel to the geometric ab— plane of a prism. The 
smallest dimension, 2c, is the sample thickness. When a 
magnetic field is oriented along the crystallographic c— 
axis (along the c— side) , the measured magnetic moment 
is denoted by Ale- Similarly, the magnetic moments mea- 
sured along the 6— side and the a— side are Mt and Ma, 
respectively. When a magnetic field is applied along the 
c— side, the Abrikosov vortices are perpendicular to the 
ab— plane and their gradients generate a supercurrent 
density, jab, given by 

- yy 20 _ 20 

'va{i-i,)--jr^ 

where V is the sample volume and Mc = ^c/V is the 
volume magnetization. This formula is written in practi- 
cal units where all lengths are in cm, magnetic moment 
is in emu and current density is in A/cm^. 

The situation is more complicated when a magnetic 
field is applied parallel to the ab— plane. Here we have 
two different current densities; one for the vortex mo- 
tion across the planes, j±, and another parallel to the 
planes, jy . By using the rectangular prism shaped sam- 
ples with substantially different a and b sides, we can 
separate these two currents. Using the Bean construc- 
tion, we have 

^''^ - 20 V Sb ) 

and, 

20 V 3a Ji,; 
Solving for the currents we obtain, 

20 bMa - aMb 

■^'-L = 1 

c — a 

and 

20 {bMa-aMbf 
~ 3a6 {Ma - Mb) (6 - a) 

where both formulas are applicable as long as aj\\ > cj± . 
If this condition is violated, similar formulas can be easily 
obtained and the applicability of a particular expression 
can be checked by examining the results. 

Fig. [13] compares the critical currents obtained from 
transport and magnetic measurements. In a broad range 
of temperatures from Tc down to 5 K, the magnetically 
determined in-plane and inter-plane critical currents are 
different by a factor of 3, slightly more than in the di- 
rect transport measurements. Fast magnetic relaxation 
in pnictide superconductors may play a substantial role 
and result in a noticeable difference between transport 
and magnetic currents (39j . Therefore, the comparative 
agreement between the two estimates is quite good. 
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DISCUSSION 



Band structure 



Comparison with other iron pnictides 



We now turn to a comparison with other iron pnic- 
tide compounds. Until now, the anisotropy has been 
characterized rehably only for 7//. Most directly this 
has been accomplished from a scaling entropy proce- 

For 
-4, 



dure in heat capacity measurements 
NdFeAs0i_2:Fj; this has given an anisotropy of -fH 
while for Bai_2:K2;Fe2As2 it is ~2.6. This is not 
very different from determinations from resistive mea- 
surements, as we discussed above. Anisotropy inside the 
superconducting state was determined from torque mea- 
surements at fields H ^ Hc2- This is sensitive to flux 
distribution in the sample (see magnetooptical imaging 
in [3^), and studies performed until now give a scatter 
of values from about 9 in SmFeAs(0,F) [53| to 1.2 in 
PrFaAs(0,F) (see Ref. [l^l for a summary of all data). 
By looking at different properties, our study reveals that 
the small anisotropy is typical for all normal and super- 
conducting properties of the iron pnictides. 

The fact that the anisotropics found in this study for 
the normal state resistivity and for the slopes of the upper 
critical fields match, implies an orbital mechanism behind 
superconductivity suppression with magnetic field. Pro- 
vided the same is the case in other superconductors in 
this family, high values of the Hc2 then would imply that 
the Fermi velocities are small, and three-dimensional por- 
tions of the Fermi surface play an important role in the 
band structure of these materials. 

It is interesting to understand whether or not a small 
anisotropy favors a higher Tc- A long history of super- 
conductivity research in low dimensional materials fol- 
lowing the original idea by Little [s^l, suggests that a 
two-dimensional electronic structure is favorable for the 
realization of a higher see for example Ref. [H^. For 
magnetically mediated superconductivity a direct link 
between the anisotropy and the Tc is suggested (56| . 
We now turn to compare the anisotropy of the upper 
critical fields in different iron arsenide materials. In 
the Co-doped version of BaFe2As2 with Te=22.5 K the 
anisotropy is about 2 from resistive 



24, m 



and 2.1 

to 2.6 from specific heat (Fig. [9]) determinations. In 
Bai_a;Ka;Fe2As2 with Tc=30K jh is in the range of 2 
to 3.5 from the resistive [iSj and equals 2.6 from spe- 
cific heat measurements [21| . In the 1111 compound Nd- 
FeAs(0,F) with Tc «50 K it is about 4 from specific heat 
[20| and TDR [l3l measurements. We may notice that 
the anisotropy is increasing in parallel with Tc- This may 
be an interesting trend that is worthy of further explo- 
ration. It actually reveals the underlying trend in the 
evolution of the band structure of these compounds. 



A key structural feature of the iron arsenides, which 
distinguishes them from the cuprates and makes a pro- 
found effect on their electronic structure, is the location 
of the As atoms above and below the layer of Fe atoms. 
We recall that in the cuprates oxygen atoms are located 
in the planes. 

Since actual location of As atoms can vary, in our cal- 
culations we used two positions. In the first one, the 
position of As atoms was determined from the minima of 
total energy. We refer to this as relaxed position below. 
The calculated position of As atom, z^s=0.341, is very 
close to that found in previous calculations, zas—OS^'^ 
[stI , but is significantly lower than the experimental value 
of 2:^s=0.355 (50I. To simulate 0.074 Co doping we re- 
placed Fe atoms by virtual atoms having Z=26.07 (vir- 
tual crystal approximation). The Fermi velocities, calcu- 
lated using the Bolz-Trap [51[ package, are very sensitive 
to the As positions. A similar trend is found in magnetic 



properties 57, [s^]. The downshift by 0.16 A increases 



both the band dispersion and the Fermi velocity along 
the z direction, see Fig. [141 The calculations with re- 
laxed ZAs give for pure BaFe2As2 and 2.7 for 
the x=0.074 doped compound studied here. Using the 
experimental zas we come to a much larger anisotropy of 
12.1 and 9.0, for pure and doped materials, respectively. 

This difference in anisotropy is affecting the most 
topology of the sheets of the Fermi surface surrounding 
r point in the Brillouin zone (Fig. I14p . Cylinders around 
X point remain warped and do not change much with 
variation of As atom positions. 

As we see, both positions of As give anisotropics not 
very different from those that have been experimentally 
determined. More precise determination of the As po- 
sition in doped materials as well as anisotropy study of 
the parent pure compound are necessary to clarify the 
situation. 



CONCLUSIONS 

Iron arsenic superconductors reveal small anisotropics 
of the electronic structure and, as a result of this, of the 
superconducting state. Since many theories consider high 
anisotropy as an important ingredient for the achieve- 
ment of high transition temperatures, this experimental 
observation puts strong constraints on the models suit- 
able for the explanation of superconductivity in these 
exciting compounds. 
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FIG. 14: Fermi surfaces of Ba(Fei-2,Coa;)2As2, a::=0.074, cal- 
culated assuming relaxed (a, top panel) and experimental (b, 
lower panel) positions of As atoms in the lattice. 
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